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ABSTRACT 

In this letter we revisit the idea of reionization feedback on dwarf galaxy formation. We show 
that primordial minihalos with v C i r < 20 km s _1 stop accreting gas after reionization, as it is 
usually assumed, but in virtue of their increasing concentration and the decreasing temperature 
of the intergalactic medium as redshift decreases below z = 3, they have a late phase of gas 
accretion and possibly star formation. We expect that pre-reionization fossils that evolved on 
the outskirts of the Milky Way or in isolation show a bimodal star formation history with 
12 Gyr old and < 10 Gyr old population of stars. Leo T fits with this scenario. Another 
prediction of the model is the possible existence of a population of gas rich minihalos that 
never formed stars. More work is needed to understand whether a subset of compact high- 
velocity clouds can be identified as such objects or whether an undiscovered population exists 
in the voids between galaxies. 
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1 INTRODUCTION 

It is widely assumed that reionization of the intergalactic medium 
(IGM) suppresses gas accretion and galaxy formation in small mass 
halos with circular velocity v C i r ^ 20 km s _1 , corresponding to 
a virial temperature T V i r , $ 20, 000 K and to dark halo masse s 
M < 10 8 - 10 9 M (e. g. . iBabul & Reesll 19921 : lEfstathiodl 19921) . 
Hereafter, we will refer to these small halos affected by reioniza- 
tion feedback as "minihalos". The reason why gas accretion and 
star formation is suppressed in minihalos is that, after reionization, 
the Jeans mass in the IGM exceeds the mass of these minihalos 
and the IGM gas is unable to condense under the influence of their 
gravitational potential . Simulations of galaxy formation have con- 
firmed this idea (e.g jBullock et al.1200 0: Gnedini |2000h. However 
a rec ently discovered dwarf galaxy - Leo T dRvan- Weber et al.l 
2008) - despite having an estimated mass of about 7 x 10 6 M 
and having properties typical of other ultra-faint dwarf spheroidal 
(dSph) galaxies, contains gas and is actively forming stars at the 
present day. This puzzling observation prompts taking a fresh look 
at the problem of gas accretion onto minihalos after reionization. 
However, the model presented in this study stands independently 
of Leo T. 

We show that the ability of dark minihalos to accrete gas from 
the IGM depends not only on the ratio of their circular velocity 
to the IGM sound speed, but also on their concentration, c. Typ- 
ically, the concentration of a halo is c V i r ~ 5 at the redshift of 



virialization but, as the halo evolves in the expanding universe, 
its concentration increases. The evolution of the halo concentra- 
tion with redshift can be understood in the con text of the theory of 
cosmological secondary infall of dark matter JBertschingeilll98 
and has been quantified us ing N-body simulations JBullock et 
l200ll : IWechsler et alj|2003) . In this letter we show that primordial 
minihalos with v C i r < 20 km s _1 stop accreting gas after reion- 
ization, as expected, but in virtue of their increasing concentra- 
tion and the decreasing temperature of the IGM at z < 3, they 
start accreting gas from the IGM at later times. As a resu lt, we ex- 
pect that pre-reionization fossils dRicotti & GnedinkOOl hereafter 
RG05) in the Local Group have a more complex star formation his- 
tory than previously envisioned. A signature of this model is a bi- 
modal star formation history with an old (~ 12 Gyr) and a younger 
( ,$ 5 — 10 Gyr, depending on the minihalo mass) population of 
stars. Another prediction of the model is the possible existence of 
minihalos containing only gas but no stars. This would revive the 
suggestion that some compact high-v elocity clouds (CHVCs) are 
extragalactic objects ( Blitz et al. 1999). However, the bulk of these 
objects may be yet undiscovered due to th eir faint 21 cm and H- 
a flu xes. The ongoing survey ALFALFA jGiovanelli etai1l2005l. 
l2007t) may be able to discover this new population of extragalac- 
tic objects. This letter is structured as follows. In §[2] and §[3] we 
describe our model for dark matter and gas accretion and in § [4] 
we apply the model to the interpretation of Leo T and observa- 
tions of CHVCs. Throughout the rest of the paper we use the fol- 
lowing cosmological parameters (h = 0.72, fi^m = 0.214, = 
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Figure 1. The overdensity profile of gas in hydrostatic equilibrium in a 
NFW dark matter halo virialized at z v i r = 0, 3, 6 and 10, (corresponding 
to halo concentrations, c = 5, 20, 35 and 50), from bottom to top respec- 
tively. The solid and dotted curves are for isothermal and adiabatic equa- 
tions of state, respectively. The /3-model provides a good fit to the isother- 
mal profile by imposing the same gas density in the core. Here, the dashed 
lines show the feeta-model renormalized by imposing that the mean gas 
overdensity in the halo is A a = 5 (see appendix). All halos have the same 
r = Tvir/Tigm = 1. The halo radii are normalized to unity. 
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2 EVOLUTION OF THE DARK MATTER 
CONCENTRATION 

N-body simulations show that the concentration of dark halos is, 
on av erage, c(z) oc (1 + z)" 1 (Bullo cket alj200ll ; IWechsler et al.l 
2002). Assuming a universal NFW density profile, the halo con- 
centration at virialization, cw, has a weak dependence on the halo 
mass (oc A'f ~ 1//6 ). If, instead, t he halo profi l e has an inner slope 
that depends on the halo s mass jRicotti|[2003l : iRicotti & Wilkinson! 
120041 : [Ricotti et af]|2007h . cw is a universal constant. Either way, 
although with some scatter and a dependence on the environment, 
the mean concentration of minihalos is c(z) = c V i r (l + z V i r )/ (1 + 
2), where c V i r ~ 5 and z v i r depends on v c i r and the cosmology. 
Using five-year WMAP cosmological parameters we have 
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where as is the variance of density perturbations at z = in a 
sphere of 8 Mpc radius, and n is a number that expresses how 
rare is the initial density perturbation that produces the mini- 
halo {i.e., the number, na, of standard deviations from the mean). 
A halo virialized from the collapse of an na perturbation has 
a typical value of the concentration c(z,M) « 50(n/3)[l — 
log (M(z vir )/10 8 M Q )" 1/5 ]/(l + z). In the rest of this letter we 
will use this equation as a rough estimate for the concentrations of 
minihalos as a function of their mass and redshift. This expression 
does not take into account the dependences of the concentration on 
the local overdensity, but, qualitatively, the model is independent 
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Figure 2. Redshift evolution of the equation of st ate (EOS) o f the low den- 
sity IGM from observations of Lyman-a forest (Ricotti et al. 2000). The 
dashed lines show the fiducial EOS used in some of our calculations. Qual- 
itatively, our results are independent of the assumed fiducial model for the 
evolution of the EOS. The data point at z = shows that the low redshift 
evolution of the IGM in regions that are not shock heated (with overdensity 
less than a few) is nearly adiabatic. 



of the assumed value of the concentration. N-body simulations also 
show that present-day isol ated dwarf halos h ave a high median con- 
centration of c ~ 35 (e.g., IColm et al.l200l) . in agreement with the 
extrapolation of the fitting formula c(z) oc (1 + Zf)/( 1 + 2), de- 
rived for more massive halos (e. g. . iBullock et alj fcoo 1 ). However, 
high-resolution simulations able to resolve the smallest minihalos 
(that are expected to have the highest concentrations) typically fo- 
cus on Milky Way type systems and only a few focus on isolated 
minihalos in the voids. 

In summary, minihalos that form at high redshift and do not 
merge into bigger halos until redshift z < 1 — 2, are likely to 
achieve a large concentration and, as a result, have a late phase 
of gas accretion from the IGM. The evolution in relative isolation 
of the minihalo also ensures that it is not tidally truncated and it 
does not evolve in the Warm-Hot intergalactic medium (WHIM), a 
phase of the IGM that has been shock heated to millions of degrees 
and believed to be the repository of most of the baryons at z — 0. 
The model has observable effects on any mini halo surviving "undi- 
gested" by a larger galaxy until 2 1 — 2 (e.g., Fig 4 shows results 
for 2-a perturbations, but similar results can be derived for l-a per- 
turbations). In particular, for relatively massive pre-ionization fos- 
sils the model predicts a bimodal star formation history. However, 
in very small mass halos significant effects are only observed in 
the rare high-er perturbations. The model that fits Leo T uses a 3-a 
perturbation. Future work will focus on estimating how many ob- 
jects like Leo T may exist in the Local Group and explore whether 
and how many minihalos smaller than Leo T could be detected as 
extragalactic compact high velocity clouds around the Milky Way. 



3 COLD GAS ACCRETION ONTO PRIMORDIAL 
MINIHALOS 

For a given value of the parameter T = (tw/Cs ffm ) 2 = 
T V i r /Ti g m, the density profile of the gas in a minihalo depends 
on its concentration and on the ability of the gas to recombine and 
cool. Figure Q] shows the gas overdensity profile for an isothermal 
(solid curves) and an adiabatic (dotted curves) gas in hydrostatic 
equilibrium in NFW halos with different values of the concentra- 
tion (the derivation of the equations for the gas density profile is 
shown in the appendix). If the gas condenses in the minihalo adia- 
batically, the overdensity in the core does not exceed 100. However, 
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Figure 3. The evolution of the gas density in the core of dark halos virializ- 
ing at z v i r = 10 with circular velocity v c i r = 18, 15, 12, 9 and 6 km s 
(from the top to the bottom). The thick curves assume an isothermal gas in 
the halo and the thin curves adiabatic. The curves are dashed if the gas can- 
not cool (tree/in > 1) an d solid otherwise {t rec /tjj < 1). The curves 
from the bottom to the top in the lower insert panel show t rec /tj{ as a 
function of redshift for the same halos as in the top panel. The EOS of the 
IGM is assume to evolve as in Fig. 2. 

if the gas condenses isothermally, the overdensity can reach values 
~ 10 7 , corresponding to gas densities of 1 — 10 cm -3 at z — 0. 

If the cooling/heating time of the gas inside the minihalo is 
longer than the Hubble time tn (that is the time scale for evolution 
of the minihalo concentration and gravitational potential), the gas 
is compressed and heated adiabatically {i.e., t coo i ~ t^eat > tn). 
As shown in Figure(2] after He II reionization at z ~ 3, due to the 
Hubble expansion, the mean temperature of the IGM decreases al- 
most adiabatically because in the low density IGM the heating time 
is theat > tn- Similarly to the IGM, during the initial phase of gas 
accretion onto a minihalo, the gas is compressed adiabatically and 
the pressure prevents the gas density in the halo core to increase 
substantially above the mean IGM density. However, if the core 
overdensity becomes larger than a critical value, the gas starts re- 
combining and when the neutral hydrogen fraction reaches a value 
xhi ~ 10~ 3 , the cooling rate increase rapidly due to Lyman-a 
cooling. At this point t coo i < tu and the gas evolution becomes 
nearly isothermal. The gas overdensity in the core of a minihalo 
can reach the large values shown in Figure[TJby the solid curves for 
an isothermal gas. 

The cooling time for a gas at T ^ 10 4 K is always t coo t < 
tree, independently of the assumed ionization fraction, tempera- 
ture and density of the gas. The cooling is due to hydrogen and 
helium recombination lines if the gas is highly ionized by UV radi- 
ation (xhi < 10 -3 ), or hydrogen and helium Lyman-a, otherwise. 
Thus, if the gas can recombine in a Hubble time, it is also able to 
cool efficiently. Figure[3]shows the evolution of the core gas density 
(top panel) and the ratio t re c/tn (bottom panel) in minihalos of cir- 
cular velocity iw for an isothermal (thick lines) and adiabatic (thin 
lines) gas. Only if t rec /tH < 1, the gas is able to collapse isother- 
mally (thick lines) and condense substantially in the minihalo core. 

The recombination time is t rec = [n ajCOre a<(T)] - , where 
a(T) is the hydrogen recombination rate. Thus, at z = we have 
trec/tn < 1 if n g , core ^ 10~ 4 cm -3 . Using Eq. JAU . we also 
have t rec /t H < 1 if0(z) > 1.5[(14/(1 + z) - 1]. Writing /3 as a 
function of v c ir we find that minihalos with v c ir > Vcir condense 



isothermally. At redshift z = 0, v% r ~ c* sm 19.5/(4.4 + c/4). 
These minihalos can be observed at 21 cm and Ha wavelengths 
and would resemble faint CHVCs. The isothermal density profile 
is well approximated by a /3-model (see Fig. [T] above) with core 
radius r core ~ 0.2r s m r V i r /25, where rw is the virial radius 
at formation. After virialization, by definition, the halo radius at 
redshift z < z vir is rhalo(z) = r vir (c/c vir ). The typical size of 
the gas core is 




although, for small values of the core overdensity, most of the gas 
extends out to r ha i « 18 kpc(M/10 8 M©) 1/3 /(l + z). The 
circular velocity, v C i r , at r core is 



v C ir(r C ore) ~ 0.66?w (r vir ) w 0.624v™ x , (4) 

where v™" x is the maximum circular velocity. 

Small minihalos that can accrete gas from the IGM may not 
be able to form stars because the gas cannot cool below T ~ 10 4 K 
if the gas has low metallicity or is metal free. These halos may 
be stellarless as CHVCs or may be pre-reionization fossils with a 
12 Gyr old stellar population, having gas but no young stars. In or- 
der to form stars the gas in the minihalo must cool below 10 4 K and 
develop a multi-phase ISM. If the gas is metal free, H2 formation 
and cooling is unlikely to be large enough to support star formation 
at z = 0. This is because, in a gas of primordial composition, the 
absence of dust grains that typically are the main catalyst for H2 
formation, require that H2 forms though a chemical reaction that 
is very slow and that involves the ion H~ as a catalyst. At z = 
the flux of the H2 photo-dissociating background (with energies 
11.3 eV < hv < 13.6 eV) that destroys very efficiently H2, is 
expected to be much larger than during the dark ages at z ~ 30. 
Positive feedback processes that may be important at high-z in in- 
creasing the form ation rate of H~ and H2 are absent at z = 
dRicotti et alj|2002h . Thus, it is unlikely that dark minihalos that 
were never able to form stars will be able to do that at z — for the 
first time. However, due to the interesting possibility of Pop III for- 
mation at z < 1 this scenario should be explored in more detail to 
asses quantitatively what is the probability of Pop III star formation 
in isolated minihalos at z = 0. 

Let's now estimate which level of metal pre-enrichment is nec- 
essary for star formation in gas rich minihalos. The cooling func- 
tion from hyperfine transitions of oxygen and carbon depends on 
the gas metallicity, Z, roughly as A23 ~ 10 -3 (Z/Zq), where 
A23 = 10~ 23 erg s _1 cm 3 . Thus, a necessary condition for star 
formation is t coo i fa (0.7 yr) T/(n 9 , core A_23) < ta, that can 
be written as n 9iCore > n* core » 0.03 cm~ 3 (Z/10~ 2 Zq)~ x . 
The left panel in Figure [4] shows n g , core and Nh in minihalos that 
evolve isothermally at T ~ 10 4 K but that do not form stars {i.e., 
candidates for extragalactic CHVCs). The horizontal lines shows 
the requirement for metal cooling and star formation assuming 
gas metallicity Z = 0.1 and 0.01 Zq. The right panel in Fig- 
ure [4] shows n 9jCore , Nh and Md yn / M gas (the dynamical mass 
to gas mass ratio) in the core of minihalos that are able to cool to 
T = 5000 K (the temperature of Leo T ISM), thus, able to form 
stars {i.e., n 9jCore > n* g core ). The symbols show the observed 
value for Leo T. For v C i r (r c ) — 7 km/s, corresponding to Leo T 
circular velocity, the model gives a core radius r c ~ 80 pc, that is 
slightly smaller than 100 pc measured for Leo T. The dark matter 
and gas mass within the core are 4 x 10 6 M and 2 x 10 5 Mq, 
respectively, in reasonable agreement with observations of Leo T. 
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Figure 4. (Left). The gas density, n g . core (solid curves), and hydrogen column density, Njj = 2r c n g . core (dashed curves), within the core, r c , of a minihalo 
as a function its circular velocity at r c at z = 0. The curves from the bottom to the top refer to 2, 3 and 4 — a perturbations. The minimum v c i r in each curve 
is determined by the condition tree/ii? < 1, necessary for cooling to T gas ~ 10 4 K. The horizontal lines with arrows show the requirement for cooling to 
temperatures below 10 4 K, necessary for initiate star formation, for gas metallicity Z = 0.1 (lower line) and 0.01 Zq (higher line). (Right). Same as in the 
left panel but for minihalos forming from a 3-<r perturbation and able to form stars due to a gas pre-enrichment to metallicity Z = 0.1 Zq and reaching a gas 
temperature T 9lC ore = 5000 K. The thick curves and thin curves refer to halos observed at redshift z = and z = 2, respectively. 



4 CONCLUSIONS & DISCUSSION 

Although Leo T contains gas, it properties as similar to those of 
other ultra-faint dSph galaxies and qu ite different from thos e of 
more massive dwarf irregular galaxies jBovill & R icotti 2001). Its 
properties are consistent with those of pre-ionization fossils but, in 
the context of this model, it is difficult to understand how it could 
hold on to its gas. 

We point out that dark minihalos with large concentration are 
able to accrete gas from the IGM more efficiently than lower con- 
centration halos of the same mass. The concentration of a halo de- 
pends on the redshift of virialization as c(z) oc (1 + 2W)/(1 + z). 
Thus, as a result of their growing concentration, early forming 
minihalos, although affected by reionization feedback, may have 
a late phase of cold gas accretion from the IGM and form stars 
(perhaps for the first time) at z < 1, This model explains all 
the observed properties of Leo T. Having a distance of 420 kpc 
from the Milky Way, Leo T is just starting to fall into the Milky 
Way halo, thus it is likely to have evolved in isolation (fulfill- 
ing a requirement of the model) and unlikely to have been tidally 
stripp ed. Recent observ ations suggest that Leo T has a bimodal 
SFH dde Jong et al]l2008h with a > 12 Gyr population and a 9 Gyr 
population (with SF continuing until few Myrs ago), in agreement 
with our expectation of a late phase of accretion from the IGM. The 
bimodality of the SFH should be most pronounced in minihalos 
with the lowest v C i r , gradually disappearing as v c %t increases to- 
ward the 20 km s _1 threshold for pre-reionization fossils (RG05). 
Leo T has a gas velocity dispersion ~ 7 km s _1 , that, according 
to our model, implies a dark halo with maximum circular veloc- 
ity of 11 km s -1 (Eq. [4j. Leo T has a core radius r c ~ 100 pc, 
M dV n /M a ~ 8 - 10 dSimon & Gehall2007l : iRvan- Weber et all 
2008) and maximum H I column density Nh = 7 x 10 cm . 
These values are all in agreement with our model as shown in § [3] 
and right panel in Fig. [4] If future observations will show that the 
dark halo of Leo T has v™£ x > 20 km s" 1 or its SFH does not 
show any decline after 12 Gyr ago, our model for Leo T will be 
disproved. 



Minihalos smaller than the one hosting Leo T can also accrete 
gas from the IGM at z < 1, but may not be able to form stars. Suffi- 
ciently small minihalos are likely to be dark and, if they evolved in 
isolation, the IGM around them would be nearly metal-free. Cos- 
mological simulations show that the fraction of primordial mini- 
halos that host a luminous galaxy decreases steeply with decreas- 
ing V dr- f(lum) ~ min[l, 50%(v c i r /W km/s) 4 ' 6 ] dRicotti et al.l 
2008). A fraction of dark minihalos may be able to accrete gas from 
the IGM but, if the gas density in their core is below the threshold 
required for star formation, these objects will resemble CHVCs. 
Additional work is needed to determine whether our model is fully 
consistent with the properties of known CHV Cs dBlitz et al]| 19991 : 
iBraun & Burtonll 19991 : kobishaw etal]|2002h . however the typical 
column densities we predict in Fig.[4jleft panel) are in agreement 
with a subset of existing observations. The typical size and size 
distribution of CHVCs may agree with the model assuming a mean 
distance of minihalos of 1 Mpc. We estimate that realistic extension 
of the diffuse gas in dark minihalos are ^5 — 10 kpc. This may 
exclude several known CHVC and/or set an upper limit for their 
distance from the Milky Way. More extended extragalactic CHVCs 
are possible for minihalo masses M > 10 s M© that are likely to 
be luminous. In addition this scenario for more massive miniha- 
los has already be en investigated by previous studies and found not 
fully satisfactory ^Sternberg et alj|2002l : iMalonev & Putmanll2003l : 
iPutman et alj|2003f) . More work in this area is also motiva ted by 
the ongoing survey ALFALFA dGiovanelli et"al]|2005l 1200% that 
may be able to discover a even fainter population of extragalactic 
CHVCs. Hence, future work shall focus on estimating the number, 
flux and size distribution of gas rich minihalos in the Local Volume 
and the local voids, producing synthetic maps for 21 cm and H-a 
wavelengths. 
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APPENDIX A: GAS DENSITY PROFILE IN NFW 
POTENTIAL 

In this appendix we derive the equations describing the density 
profile of intergalactic gas inside a dark minihalo. Although, the 
IGM density, temperature and dark halo potential evolve with red- 
shift, the assumption of quasi-hydrostatic equilibrium is a good 
approximation within the halo radius. Expressing the halo radius 
as r ha i = v cir t H /(4:A dm ) 1/2 , where A dm ~ 178 is the halo 
overdensity and tn is the Hubble time, we find that t cros < tn 
(where t cr0 s = rhaio/cl 3m is the halo crossing time) if T = 
{v clr /c l s 9m ) 2 < (4A dm ) ~ 800. In this letter we are interested 
in halos with V 1, thus t croa < tn and we can safely assume 
quasi-hydrostatic equilibrium in order to calculate the gas density 
profile. We assume a NFW density profile, p oc l/[cx(l + cx) 2 \, 
with x = r/rhaio, halo concentration c = rhaio/r 3 , where r„ the 
core radius of the dark matter profile. The overdensity profile of a 
gas with equation of state P — Kp 1 , in hydrostatic equilibrium in 
a NFW halo of mass M, concentration c and circular velocity, v c i r , 
is 

1 + 5 g (x) = f- = [1 + ( 7 - 1) In (1 + 5i s °)] ^ (Al) 

Pb 

where pt is the mean IGM gas density and 1 + 5 g so (x) — (1 + 
g 

cx) ™ is the overdensity for 7 = 1 (isothermal equation of state). 
In Eq. lAll we have normalized the density profile so that S g — 
at x — > 00. The parameter /3 = Tc/f(c) ~ T(4.4 + 0.25c) de- 
termines the core overdensity (at x = 0). The gas density profile 
has a core with overdensity 1 + 5 g , core = [1 -+- (7 — l)/?] 1 ^ 7-1 ^ , 
that for 7 = 1 is 1 + 5 g s ° COTe = e 8 . When S g a ° ore > 10 6 - 10 7 , 
the normalization of the density profile becomes inaccurate be- 
cause the mass inside the halo exceeds the cosmic value within 
the turnaround radius. In this regime the density profile must be 
normalized by imposing that the total gas mass within the halo 
equals a given fraction of the dark halo mass (this is necessary be- 
cause the calculation of the gas density profile does not take into 
account the Hubble flow and would overestimate the gas accre- 
tion from the IGM in this regime). The /3-profile 1 + 6g(x) = 

5 C [l + (5cx) 2 ] 36 2 , where b — 0.2/3/3, provides a good fit to 
the isothermal density profile and has the advantage of being rela- 
tively easy to normalize to a constant mass M gas = 47r/3A 9 p6r^, 
where A g ~ 5 is the mean gas overdensity inside the halo. The 
core density for the /3-profile with the aforementioned normaliza- 
tion is 5 C = A 9 (l - 6)ai/(a< 1_i,) - b), with 01 = (c/0.2) 3 . We 
will use this expression in our calculations of the core density for 
7 = 1. 
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